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There are two aspects to the statistical physics of systems with quenched disorder. One is 
the static aspect, namely to understand the properties of the Boltzmann Gibbs distribution, 
. the other is the dynamical problem. In these lectures I would like to concentrate on some 
^ ■ recent developments about the statics. The interested reader can find some brief description 
of the recent works on dynamics in [|TJ. A system with quenched disorder contains two 
types of variables, q and t, which equilibrate on very different time scales. The thermalized 
variables t are assumed to reach thermal equilibrium (this is an assumption of the static 
problem, it is in general not true for real materials, and then the dynamical approach is 
S j needed at low temperatures). The quenehed variables , are given a priori and do not 

thermalize. In the problems we shall study there are many quenched variables, which means 

Ch ! 

a number which diverges in the thermodynamic limit where the number of quenched variables 

O 



X 



diverges. One sample of the disordered system corresponds to one set of values of q. We 

would need an infinite amount of information to just describe one sample. So in practice 

what we know is the probability distribution, "Pflll, of q. 
03 ' 

In practice there are many systems in this category which have been studied in recent 
years. A canonical example is the spin glass problem [0|§. The techniques and ideas from 
the spin glass mean field theory have been applied to many other systems, like directed 
polymers (see for instance 0-0). One could also mention the Random Field Ising Model 
(RFIM) fUJ, the protein folding problem |To| , pT| , and several problems outside of the usual 



physics ones, i.e. where the 'energy' function is not a real physical energy, like neural 
networks (see for instance jl2fl), or optimization problems (see for instance 0). 

Thermal equilibrium means that the probability of a given configuration of thermalized 
variables is given in terms of the energy H[q,t] by the Boltzmann weight: 

^H^jX>-™ (i) 



1 



where Z[q] is the partition function. So there is one such probability distribution for each 
sample. We face a difficult problem which is to characterize the set of the P q [t\. Fortunately 
it turns out that some quantities are self averaging. This means that they become sample 
independent (i.e. independent of q) in the thermodynamic limit, a typical central limit 
theorem. This is the case for thermodynamic quantities, and it is easy to derive this result for 
finite dimensional systems with finite range and bounded interactions. The properties which 
are not self averaging may also be interesting (for instance the nature of the ground state 
of the travelling salesman problem, or the folding of a protein), but this is an algorithmic 
problem. Here we shall keep to analytic studies of self averaging quantities and properties. 

One standard approach of statistical physics uses mean field methods to work out the 
phase diagram of the system. This is efficient away from the second order phase transitions. 
Around such transitions the renormalisation group can be used. It has been gradually 
understood that in many disordered systems this very general and powerful approach fails. 
The fundamental reason for this failure is that the free energy functional possesses many 
secondary minima (metastable states) which can be quite different from the minimum, and 
are not taken into account in the standard approach. A striking example is provided by the 
RFIM. It has been shown that, order by order in perturbation theory the critical exponents 
of a d dimensional Ising system with a quenched random field are equal to those of the pure 
system in dimension d — 2 This is one of the few results in field theory which is known 
to hold to all orders in perturbation, and a beautiful derivation can be obtained through 



supersymmetry |TJ]]. However this result differs from the simple domain wall argument of 
Imry and Ma \TE\, and it would predict a lower critical dimension equal to 1, while there 



are by now exact results which show that this lower critical dimension is equal to two |16 



|L7| . It has been realized long ago that the failure of the perturbative approach is related to 
the existence of many metastable states in this problem |18| . 



What one needs is a systematic and non perturbative method which is able to handle 
the field theory of disordered systems with many metastable states. Such a method can be 
obtained with the use of replicas, to handle randomness, together with a gaussian variational 
method. It turns out that the possibility of breaking replica symmetry allows to handle the 
metastable states. Such an approach was used on the problem of heteropolymers or proteins 
T~0j , p7TJ1 , but this is a complicated problem which requires additional assumptions |19 |. It 



has been developed as a systematic field theoretic method to handle manifolds in random 
media in and used since on several problems. Before turning to this case, let us 

just comment on the general meaning of replica symmetry breaking (rsb). The replica 
method allows to compute extensive thermodynamic quantities like the free energy. Using 
the property of self averageness, we need to compute the average of the logarithm of the 
partition function, which is written as: 
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ln(Z) = V(1I) ln(Z(H)) = lim ^ °° (2) 

The n'th power of the partition function can in turn be written as the partition function 
for n replicas of the original sample (with the same disorder). Disorder averaging leads to a 
pure problem (no disorder), with an effective attraction between the replicas: 



z( ff )» = 5>(n)z(n)\ = 

M 

= £P(H) ]T e W -i3(n[U,U ao } + ... + H[U,U\})= £ exp ...,U\ 

[5] Uoo,...,U\ Uoc,...,U\ 

(3) 

The replica Hamiltonian H e f/[tx, ■■■,t n ] is obviously invariant under permutations of the n 
replicas. The problem is to know whether this symmetry is preserved or spontaneously 
broken in the thermodynamic limit. As first found in spin glass theory (for a review see ||), 
such a spontaneous symmetry breaking can occur, and is related to the appearance of a spin 
glass phase. Roughly speaking, if the system can have generically many metastable states, 
each replica can condense in one or the other such state, leading to this rsb phenomenon. 
The problem is in the term 'generically': In order for Z(q) n to be typical of the generic 
behaviour, we need that n go to zero, and therefore this rsb requires understanding the 
properties of the permutation group when n — > 0. This was achieved by Parisi in 1979 and 
is reviewed in J3[. A more rigorous definition of rsb can be found using two real replicas 



(i.e. n — 2), coupled through an extensive energy term of strength e |20|. We consider two 
replicas with the same disorder governed by the Hamiltonian: 

H e [q, tuh] = H[q, h] + H[q, t 2 ) - e5[h - t 2 ) (4) 

(Here the last term is written in a symbolic way. It will in fact depend on the problem at 
hand. What I mean is an extensive term which is the integral of a local interaction between 
the two replicas, and the interaction is attractive if e > and repulsive otherwise). Calling 
F(e) the quenched average of the free energy for this two replica system, the onset of rsb 
is signaled by a non analyticity of F(e) at e = 0. The idea is the following: the difficulty 
in random systems is that we do not know what is the conjugate field which selects one 
given state (such a knowledge depends on the sample and requires an infinite amount of 
information). In the above procedure each replica plays the role of the conjugate field for 
the other one. Practical implementations of this idea can be found in for spin glasses 
and in [0] for directed polymers. 
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II. MANIFOLDS IN RANDOM MEDIA 



Consider a D dimensional manifold described in the solid on solid approximation by a N 
component vector field u(x) where x is the D dimensional vector of internal coordinates. The 
Hamiltonian describing the system is the sum of a rigidity term and an external potential: 

If D ( duo \ 2 r 
h[u] = -J dx Xj y-^r) + J dx V{x,uj{x)). (5) 

The external potential is random and gaussian. Its correlations are given by [p3fl : 



V(x, cu)V(x', uj>) = -5 {D \x -x')N f r- ^ j . (6) 

In most of the cases the relevant part of the correlation function of the potential is its 
asymptotic behaviour for large transverse distances, which we suppose to be described by a 
power law: 

,2\ (1-7) 




g uj 



2(1-7) \N t 



(7) 



One rather general class of problems (for a review see for instance |2l|j24l] ) consists in un 



derstanding the transverse fluctuations of the manifold, governed by the exponent ( defined 
by: 



< [uj(x + 1)-uj(x)}*> ~ l>>a f(. (8) 

The general model is mainly described by three parameters: D is the internal dimension 
of the manifold, N = d — D is its codimension, and 7 characterizes the large distance 
correlations of the potential. For N — 1 the manifold is an interface, for D — 1 it is 
a directed polymer, for N = D it describes the elastic deformations of a D dimensional 
crystal with impurities, and for D = 3, iV = 2 it describes the elastic deformations of a 
vortex lattice in certain pinning regimes of type two superconductors f26|. 



III. THERMAL FLUCTUATIONS WITHOUT DISORDER 

If V — 0, the fluctuations are purely thermal. The propagator is 1/k 2 , leading to ( = 
if D > 2 (the manifold is flat), and ( = (2 — D)/2 if D < 2 (the manifold is rough), and all 
these results are independent from N. 
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IV. RANDOM FORCES 



Another simple case is that of random forces, where the pinning term of the Hamiltonian 
is — / dx f(x)uj(x), with gaussian random forces of strength F. This may be relevant for the 
deformation of a cristal lattice on small enough length scales so that the displacement of an 
atom is much smaller than the correlation length of the potential. The problem is simple 
because there are no metastable states. One gets immediately 

< >= f ~§- (9) 

so that: 



< uj a (k) >< uP(k') > = 5 a/3 F 2 5(k + k') — . (10) 

K 

Clearly the l/k A term leads to a ( exponent which is zero above D — 4, and equal to 
(4 — D)/2 in 2 < D < 4 (In D < 2 the lateral fluctuations between two points at distance x 
depend on the total size of the manifold). The thermal fluctuations are much weaker since 
the connected correlation function is: 

T 



<(u«(k)-<Lj«(k)>y> = — . (ii) 



The lesson from this computation is that the disorder fluctuations are dominant, and they 
lead to a typical "dimensional reduction" result, namely the fact that the value of ( for the 
disordered system in D dimension is equal to the value of ( for the thermal fluctuations of 
the pure system in dimension D — 2. This is exact for random forces. 

V. RANDOM POTENTIAL: VARIATIONAL APPROACH 

This is a complicated problem with many metastable states. A simple argument is that 
of Imry and Ma ||15|1 . Let me just do it in the case of an interface (N = 1) in the RFIM, which 
corresponds to 7 = 1/2. A bump of size u ~ w on a length scale L costs an elastic energy 
L D ~ 2 w 2 , while the typical gain from the pinning energy is of order V L D w. Optimizing w 
leads to w ~ £( 4 - D )/ 3 ; so that £ = (4 — D)/3. Generalizing this argument to arbitrary 
correlations of the potential, we find: 

C = WTY,' (12) 

This is a reasonable guess for the exponent (, but it is difficult to improve it and it has no 
reason to be exact (a well known counter example is the case N = D = 1,7 = 3/2, where it 



is known |]33| , |32|| that ( = 2/3, while the above formula leads to ( = 3/5). 
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What one would like is a general field theoretic approach similar to that of the previous 
sections to handle this problem. A first step is to use perturbation theory. Although this is 
not strictly necessary, it is convenient to introduce replicas. The effective replica Hamiltonian 
after disorder averaging is: 

1 r « *> (duo a ^ 2 



H„ — — 
2 



/^EE3 +f/<k EM*)) 2 (is) 

J a =l fj,=l \ UX ft/ Z J a=l 

+f/*s*/( Mj) ~'» (J)1 ' y < 14 > 

a,b \ I 

The usual perturbative treatement of the Hamiltonian H n consists in expanding / in powers 
of the quadratic part gives the free propagator F ab : 

p m tab g) n ^ 

abl J k 2 + n2f3f'(0) k 2 (k 2 + n2f3f'(0)Y 1 J 

For n — ► , the presence of the 1/k 4 term immediately leads to £ = (4 — D)/2, as in the 
random force case. This result holds at higher orders in perturbation. We face a typical 
situation where perturbation theory, which is unable to handle the metastable states, gives 
qualitatively wrong results. 



The method developed in |24J uses a variational method which evaluates the best 
quadratic Hamiltonian to compute the quenched free energy with the replica method. I 
shall describe only the main ideas here, leaving all the (sometimes complicated) technical 



steps aside. The reader is referred to pi for a full presentation. The case D = 0, studied in 



[TJ is a somewhat simpler exercise which is also of interest. The most general quadratic 
Hamiltonian is: 

H v = - J dk G~\(k) uj a (x)uj b (x) (16) 
Optimizing the variational free energy leads to a simple gap equation for the self energy a a t,: 

G~l{k) = k 2 5 ab - a ab (17) 

a ab = 2(3 f ^ J dk[G aa (k) + G bb (k) - 2G ab (k)] ] j , a^b (18) 

<?aa = ~ E a ab- ( 19 ) 

It is easy to seek a replica symmetric solution to these equations, i.e. to assume that 
&ab — One immediately finds again the famous 1/fc 4 term in the propagator, leading 
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to ( = (4 — D)/2. This result is wrong mathematically because this saddle point is not a 
local minimum of the variational free energy. It is also wrong physically The reason is the 
following: The most general replica symmetric quadratic Hamiltonian describes a situation 
in which: 

du a ex P {-^ J dk ffc 2 $> a (fcK(-A:) -cr5> a (fcV 6 (-fc)j j . (20) 



Disentangling the square, one gets: 



d[f)exv L±- I dk\f(k)\< 



duexp ( J dk (k 2 \u(k)\ 2 - Lo(k)f{-k)) 



(21) 



In this form it is clear that the situation we describe is nothing but a random force problem, 
with a quenched random force with a distribution exp (— ^ J dk\f(k)\ 2 ^j. The replica sym- 
metric quadratic Hamiltonian is bound to describe the random force problem, a problem 
without metastable states. To get better results one needs to break replica symmetry. 

The rsb solution has been worked out using the Ansatz developed by Parisi for the spin 
glass problem (see the review in [[!]). This Ansatz describes hierarchical n * n matrices 
in the n — > limit by continuous functions on the interval u G [0,1]. Therefore the self 
energy becomes a function of this internal parameter u. We get a spectrum of masses, which 
results in a non trivial behaviour of the correlation function. The wandering exponent turns 
out to be identical to the one derived above through the Imry Ma argument. Therefore it 
provides a microscopic "derivation" of this result. I shall not reproduce here these technical 
computations, but rather sketch in the next section the physical interpretation of the result. 
Let me point out a few properties of the gaussian rsb variational method. It is a very versatile 
method which can be used on many problems. It is basically a Hartree approximation 
which becomes exact in the large N limit. A systematic expansion around iV = oo is thus 
in principle possible, although technically difficult. Such an expansion was carried to first 
order on the RFIM, because the leading Hartree approximation had no rsb (see ||). In the 
case of the random manifold, some attempts have been made [p9[|34"| , but the full expansion 



has not yet been worked out. Only recently has the replica technical apparutus needed for 
this expansion been developed |35| . 



VI. PHYSICAL INTERPRETATION OF THE SOLUTION 

It is very instructive to work out the physical content of the rsb solution, in the same 
spirit as we did for the rs solution before. The idea is to consider seriously the quadratic 
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Hamiltonian with rsb, and deduce what kind of physical situation it describes. This is not 



an easy task. Again I will mention only the results, referring the reader to for the 
derivations. Let us concentrate on one degree of freedom uj = uj(x) — o>(0). 

We have already seen before that in the rs case uj feels a random force. It means that 
for a given sample there is one value of the force /, and the probability distribution of uj is 
c* exp(-a(^ - ff). 

The next stage of approximation uses a single breaking step in Parisi's hierarchical con- 
struction. Its physical interpretation is as follows: For each sample one generates a set of 
favoured values of u, called uj a , together with a set of weights W a . These variables uj a and 
W a are quenched random variables. The favoured values uj a are independent variables, with 
a distribution 

2^ 2 Ql 



= 7^72 ex P -S 1 ) (22) 



As for the weights W a , they are derived from some "free energy" variables f a through: 



e -Pf a 

The f a are independent random variables with an exponential distribution such that the 
average number of states with free energy less than / (i.e. weight W a greater than e - ^) is: 

Af(f) ~ e?f, (24) 

For a given sample, that is given the variables uj a and W a , the Boltzmann probability for uj 
is given by: 

P(u;) = <? £ W a exp (25) 

This distribution of uj is characterized by three numbers qo,qi,p, which are the variational 
parameters which are determined by the variational method at this one step rsb stage. The 
above formula (|25[) is somewhat surprising. It parametrizes the Boltzmann distribution of 
the degree of freedom uj by a weighted sum of gaussians. Such a representation clearly allows 
for the existence of several metastable states contributing simultaneously to P(uj), which is 
nice. The surprising point is that such a situation is possible using a gaussian approximation. 
The power of the method comes form the combined used of the gaussian variational method 
together with the breaking of replica symmetry. An elementary example would be to look 
at a particle in a symmetric double well potential, and approximate its Boltzmann weight 
by a gaussian. At low enough temperature the best gaussian is shifted from the origin 
and breaks the inversion symmetry in the problem. To restore it one approximates the 
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Boltzmann distribution of the particle by the sum of the two symmetric shifted gaussians, 
and one finds a pretty good approximation to the exact distribution at low temperatures. 
The situation is exactly analogous here, but the symmetry at hand is the replica permutation 
symmetry. 

The full construction consists in iterating the above procedure. For two steps of rsb one 
must divide each gaussian in (p5|) into a weighted sum of subgaussians, introducing two new 
parameters (one 'q' parameter for the width of the subgaussians, one 'p' parameter for the 
distribution of their relative weights). The construction is then iterated an infinite number 
of times, adding each time two new variational parameters. Clearly in this way we can 
generate very complicated probability distributions for u>. 

To summarize, the formalism of broken replica symmetry provides a large space for the 
probability distribution of P(uj) (when one changes sample). The probability distributions 
so generated have the following properties: 

a) They depend on many parameters so that we have a large variety of choices. 

b) If we want, we can construct these probabilities in such a way that the system is 
scaling invariant at large distances, characterized for instance by a very non trivial scaling 
exponent (. 

c) The expectation values can be computed explicitly so that this solution can be taken 
as the starting point of a perturbative expansion. 

d) Last, but not least, the gaussian rsb Ansatz becomes exact when the dimension of the 
space goes to infinity. 

An interesting example of application of this very general scheme to the vortices in 



superconductors can be found in [26 
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